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Abstract 

[_L] ■ Pulsars are believed to be magnetized neutron stars. Their surface magnetic 

field ranges from 10 8 to 10 12 G. On the other hand, the magnetars have 
surface magnetic field 10 14 — 10 15 G. It is believed that at the center, the 
magnetic field may be higher than that at the surface. We study the effect of 
the magnetic field on the neutron star matter and hence on the mass-radius 
relation of neutron stars. We model the nuclear matter with the relativistic 
mean field approach considering the possibility of appearance of hyperons 
at higher density. We find that the effect of magnetic field on the matter 
of neutron stars and hence on the mass-radius relation is important when 

i/-) ' central magnetic field > 10 17 G. Moreover, if the central field is of the order 

of 10 19 G, then the matter becomes unstable which limits the maximum 

^+ . magnetic field at the center of magnetars. 
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1. Introduction 

Soon after the discovery of radio pulsar [l|, the theoretical proposition of 
neutron star (NS) [2j drew much attention. At the core of a NS, the matter, 
composed of neutrons, protons, electrons, and sometimes muons, is highly 
degenerate with density in the range of 10 14 — 10 15 gm/cc. Naturally, the 
constituent particles therein interact via nuclear forces, forming a non-ideal 
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nuclear fluid. However, the nature of nuclear forces at this high density and 
the internal composition, especially the composition of matter in the inner 
most core of the NS, are not well understood yet. There are many theoretical 
models of nuclear matter at high density describing different equations of 
state for the matter. Based on different equations of state, different mass- 
radius relations of NS are obtained. Consequently, only way to constrain 
the equations of state is comparing the theoretical results with observed 
properties of NS. 

With the detection of pulsars, it was also observed that the pulse pe- 
riod slowly increases in a very regular manner. Consequently, pulsars were 
successfully modeled as magnetized rotating NSs [3, |4j, [5[]. This model suc- 
cessfully describes the spin-down power of the Crab pulsar and other pulsars. 
This description requires a surface magnetic field 10 11 — 10 13 G for the pul- 
sars detected at the beginning. Later on, more observations broadened this 
range of the magnetic field with the discovery of millisecond pulsars [6j . This 
class of pulsars has a much lower surface magnetic field of 10 8 — 10 10 G [7J. 
In the recent past pulsars with a very high rate of change of period (period 
derivatives) have been observed |8j,[9|], implying a high surface field up to 10 14 
G. Another class of pulsars has been observed which is the accreting X-ray 
pulsar 10(| . For this class, the surface magnetic field varies from 10 8 to 10 13 



G. 

Anomalous X-ray pulsars (AXPs) were discovered in early 80s, which 
show periods in a relatively narrow range of 5 — 11 s and have large period 
derivatives [llj. In early days, they were believed to be unusual accreting 
NSs, whereas, later surveys detected no companions for these candidates. 
The large value of period derivatives implies a surface magnetic field as high 
as 10 14 — 10 15 G. This helps in explaining the X-ray luminosities, powered by 
magnetic field decay, without any accretion. However, until the discovery of 
soft 7-ray repeaters (SGRs), such high values of surface magnetic field were 
not acceptable. SGRs show brief intense bursts of soft 7-rays. From the 
peak luminosities of the bursts, SGRs are identified as NSs with very strong 



surface magnetic field (~ 10 14 - 10 15 G) [l2|,[l3(. Now both AXP and SGR 



are believed to be strongly magnetized NSs with different variations, which 



are popularly known as magnetars [11|, [14], [15|, [16 . 

The existence of a magnetar motivates to study the effects of strong 
magnetic field on NS properties. A strong magnetic field affects, the structure 



of a NS through its influence on the underlying metric [17|, [18( and equation 



of state (EoS) through the Landau quantization of charged particles and 



then the interaction of magnetic moments of charged particles (and even the 
anomalous magnetic moment (AMM) of neutral particles) with the magnetic 
field. In the present paper, we plan to study the effect of magnetic field on the 
EoS of highly dense matter. Generally, the effect of magnetic field on the EoS 
of matter is significant when the magnetic field B > 10 18 G. For the nuclear 
matter with a n-p-e system, the effect of magnetic field was studied by several 



authors |19| . |20L |2lL l22l |23| . However, the composition of the core of a NS is 
very uncertain as the density reaches 10 — 15 times the normal matter density. 
At this high density, the quarks may be deconfined, making quark matter, 
or the hyperons may appear, making hyperonic matter. Another hypothesis 
is that the whole star is made of deconfined quark matter, which is known 
as strange star (SS). The effect of magnetic field on quark matter using the 



MIT bag model has been studied in different literature [24j, |25|, |26| . There 



are other models of quar k matter with phenomenological density dependent 
quark masses [27|, |28|, |29|, |30( . Broderick et al. [31] studied the effect of strong 
magnetic field on hyperonic matter. In their study, the field strength does 
not depend on density. However, in realistic situation, it is believed that 
the field strength is higher at core than that at surface of a NS. Therefore, 
with radius, and hence with density, the field strength should vary. In the 
present work, we plan to study the influence of the strong magnetic field 
on the hyperonic matter and hence on the mass-radius relation of NSs with 
hyperonic matter in the core with a density dependent magnetic field. 

The paper is organized as follows. In the next section we describe the 
model of hadronic matter at high densities in the presence of magnetic field. 
Subsequently, we discuss the numerical results in section |3j Finally, we sum- 
marize the results in the last section. 



2. Model 

Among many theoretical models of dense nuclear matter, the relativis- 



tic mean field theory of nuclear matter, also known as Walecka model [32 



successfull y d escribes the nuclear ground state properties and elastic scat- 
tering 



32, 33, 34, 35 



In this model, the nucleons in matter interact among 
each other through three meson fields: the isoscalar-scalar meson a, isoscalar- 
vector meson u, and isovector-vector meson p. However, the original Walecka 
model 32] predicts very high incompressibility of nuclear matter. One of the 
ways to solve this problem is by introducing cubic and quartic self interaction 
terms for the scalar field into the Lagrangian 36|. With this prescription, 



the EoS of nuclear matter and the properties of NSs were studied extensively 



37, 38, 39, 40, 41 



In the present paper, we employ this idea to construct the model of dense 
matter. In addition, we consider the possibility of appearance of hyperons 
(A, E~, S°, S + , H - , S°) and muons (/i~) at higher density. The Lagrangian 
density of matter under consideration is then given by J3j 
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Here, B indicates all baryons present in the matter. The scalar self interac- 



tion term 



is 
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with the choice of gauge A^ = (0, — yB, 0, 0), where B being the magnitude of 
magnetic field and eQ the charge of the particle, e the positive unit of charge. 
For this particular gauge choice B = Bz. This does not mean any preferred 
choice of direction, but only specifies the coordinate axes. ^>_b, -0$, <r, u; and p 
are fields of baryons, leptons, a-mesons, w-mesons and p-mesons, with masses 
m B , mi, m a , m u and m p respectively, g aB , g wB and g pB are coupling constants 
for interactions of a, u and p mesons respectively with the baryon B. In the 
mean field approximation approach, the baryons acquire the effective masses 
as 

m B = m B - gaBcr, (6) 



where a is given by its ground state expectation value 
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where p B is the momentum and p F the Fermi momentum of the baryon B. 
In the presence of magnetic field, the motion of the charged particles is 
Landau quantized in the perpendicular direction to the magnetic field. With 
the given gauge choice for B described above, the positive energy solution 
for any charged particle is given by 



^ = g-i^t-p.X-Jrfl^^^p^ 
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where p x , p z are the components of momentum along x and z directions 
respectively, E n the energy at the nth Landau level and U s the suitable 
spinor. For positively charged baryons, up and down spinors are given by 
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and for negatively charged baryons, 
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The spinors of negatively charged leptons are same as that of negatively 
charged baryons except m* B replaced by mi which is the mass of the leptons. 
Here 
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with H n (C,) being the Hermite polynomial. The momentum in the x-y plane 
is quantized and hence the energy in the nth Landau level is given by 



E n = ^ V l + m 2 + 2ne\Q\B. 



(16) 



For positively charged particles, e.g. proton, n = 0,1,2... for the spin up 
states and n = 1,2,3... for the spin down states. However, for negatively 
charged particles, e.g. electron, n — 0, 1,2, ... for the spin down states and 
n = 1,2,3... for the spin up states. The number density of neutral baryons 
is given by 
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and the number densities of charged baryons and leptons are given by 
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where n max is given by 
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The scalar density of neutral baryons is 



n ( B ) = mf hj^ _ yf+^ \ (20) 

S,7V 7T 2 V2mf 2 ■ ' ' 



B ~ m B 



while for charged baryons it is 



b) e\Q\B * v^* , PB{n)+fi* B 

ic r " -^T m B 2^ ( 2 ~ ^°) ln / = = ' ( 21 ) 

n=o \/ m B + 2ne|(5|B 



Z7T 



with 



A*b = yJPF )2 +™*B andp B (n) = y^ 2 - 2ne|Q|£. (22) 

The chemical potential of baryons is hb = /•*# + ^Pg^B + p-j^ y^s where 
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which are the ground state expectation values of uj and p% fields, and 1% is 
the isospin projection of the baryon B. The total baryonic number density 
of the matter is rib = ^2 B Ub- The matter is neutral in charge, which implies 
n p +n^+ = n s - +ns- +n e +n^. In addition, the matter is in beta-equilibrium, 
which gives the condition /i n — /i p + \i e . 

The total energy density of matter is then 



J- o o T w \ -L on2 1 



-m 2 y + U(a) + -mluj + -m 2 pP l 

El / Cjv) *3 (AT) * 2 * * 4 i 

—2 2p F 'n N - p F J m N fi N - m N In 



iV 



Pf + A* AT 



m^ 



e\Q\B 
(2vr) 2 



EE( 2 - m 



C n=0 

x [PcH/*c+ ( m c + 2ne|Q|£)ln 
e|Q|B 



Pc(n)+/^c 



yVc + 2ne|Q|i3) 



(27r)2 EE( 2 -^ 



x [pi{n)fjn+ (m 2 + 2ne\Q\B)ln 
B 2 



Pi{n) + Hi 



^(mf + 2ne\Q\B) 



(24) 



where N is charge neutral baryons, C charged baryons and pi(n) = \J p F — 2ne\Q\B. 
The total pressure is then 

P = y] P'BUb + ^ pint - e. (25) 

B I 

The relation between energy density and pressure describes the EoS of the 
matter. 

Now we employ this EoS to get the mass-radius relation of a NS with 
hyperon in its core. We consider the stationary and axisymmetric metric, 
which can be written as 42 
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where u,j, a and w are functions of r and 9 only. Next, we construct and 
solve the equations of hydrostatic equilibrium following the previous authors 
42. 



3. Results and discussion 

Nucleon-meson coupling constants are chosen in such a way that the 
nuclear matter properties can be reproduced as the binding energy E/B = 
— 16.3 MeV, the saturation density n Q = 0.153 fm -3 , the asymmetry energy 
coefficient a aiw =- 32.5 MeV and the incompressibility K = 240 MeV, and 



asy 



are taken from previous work 37]. The hyperon-o; coupling constants are 



determined from SU(6) symmetry of the quark model |43l . |44|, |45[. The 
hyperon-o" coupling constants, however, are determined from the potential 
depth of hyperons in normal nuclear matter given by 

U Y = -g^ycy + gwYUo, (27) 

where Y corresponds to hyperon. We take the potential depth for A- hyp eron 



Ua = —30 MeV, as obtained from the analysis of A-hypernuclei |43|, |46 



From various experimental data [46|, |47J for S-hypernuclei, we obtain the po 



tential depth for S-hyperon Us = —18 MeV. Recent S-hypernuclei data in- 



dicate that the corresponding potential depth to be repulsive [48J and we take 
Uy, = 30 MeV . Note importantly that earlier work [3l| assumed hyperon- 
meson coupling constants to be fixed fractions of respective nucleon-meson 
coupling constants and for all the hyperons these fractions are same for a 
given meson. Hence, unlike our model, in that work the first appeared hy- 
peron is E + , which does not appear at all (see Tabled]) in our model. 

In a realistic situation, the magnetic field at the core may be higher than 
that at the surface and gradually may decrease towards the surface. For the 
magnetars, the maximum surface magnetic field is 10 15 G. For this case, the 
core magnetic field may be of the order of 10 17 — 10 18 G. Without any proper 
knowl edg e of magnetic field configuration inside a NS, following previous 



work |49| here we adopt a magnetic field profile 



B{n h /n Q )=B a +BAl-e- fi \%) |, (28) 

where (3 and 7 are two parameters determining the magnetic field profile with 
given B s and B c , and n& is the total baryon number density. We consider 
different values of field at the center, whereas surface field strength is taken 
to be B s = 10 15 G. We notice that the different values of B s do not affect the 
EoS considerably, and thus we do not show the results with other B s . In the 
above parametrization, the magnetic field strength depends on the baryon 
number density. However, at each density the field is uniform and constant. 

3.1. Equation of states 

We compute the EoS self-consistently with the condition of beta-equilibrium 
and charge neutrality. It is found that, as described by the Figures given be- 
low, the effect of magnetic field is important for B c > 10 17 G. However, for 
B c < 10 18 G, the effect is significant only when the field reaches B c at a very 
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Table 1: Threshold densities for muons and hyperons to appear in units of hq. 



low density, away from the center, and remains almost constant up to center. 
Therefore, we restrict our study with B c ~ 10 18 G. At lower densities, the 
matter is composed of only neutrons, protons and electrons. Hence, at the 
low density regime, the particles which are affected by the magnetic field are 
electrons and protons. Since the electrons are highly relativistic, electron 
Fermi momentum is very large compared to electron mass. Therefore, the 
number of occupied Landau levels by electrons is very large, eventhough the 
field strength under consideration is larger than the critical field strength 
of electron by several orders. On the other hand, the field strength under 
consideration is very less than the critical field strength of protons. Conse- 
quently, the number of occupied Landau levels by protons is also large. As 
density increases, the heavier particles appear gradually. In addition, the 
magnetic field increases with the increase of density. As a result, the number 
of occupied Landau levels gradually decreases for every species. The thresh- 
old densities for muons and other hyperons to appear for B c = 10 18 G are 
given in Table HJ The threshold densities for various species to appear do 
not differ from their respective values when the magnetic field is absent. 

In Fig. [1] we show the EoSs with and without magnetic field when the 
magnetic field profile corresponds to j5 — 0.1 and 7 = 1 [see Eq. (|2"8]) ]. We 
observe that the EoS becomes softer with the increase of E c . Here we should 
mention that, when the field strength is high enough, the field energy and 
the field pressure are not negligible. In calculating the EoS, we thus add this 
contribution too which is necessary to construct the structure of NSs. 



3.2. Mass-radius relations 

With the above EoSs, the structure of NSs with and without rotation is 
constructed. In this preliminary study of NS structure, we do not include 
anisotropic solutions based on the general relativity which is important in 



10 



300 




n b /n 

Figure 1: Variation of P as a function of normalized baryon number density. The solid 
curve is for without magnetic field. The dashed, dotted and dot-dashed curves correspond 
to B c = 10 18 , 2 x 10 18 and 3 x 10 18 G respectively. 



strong magnetic field. Without any rotation, the Einstein's equation leads 
to the Tolman-Oppenheimer-Volkoff (TOV) equation 



dP 
dr 



[e(r) + P(r)] [M(r) + 4vrr 3 P(r)] 



r[r-2M(r)} ^ ' 

where r is the distance from the center of the star and M(r) the mass enclosed 
within r. The radius of a star (R) is the distance from the center at which 
the pressure becomes zero. For different central densities, stars with different 
masses and radii are formed. In Fig. [2] the mass-radius relations for static 
stars which correspond to EoSs described in Fig. [I] are shown. Since the 
effect of magnetic field softens the EoS, in the presence of magnetic field 
the maximum mass is less than that for without magnetic field. Stars with 
rotation, as expected, have larger radii than their static counterpart. In Figs. 
|3] and H] we show the corresponding mass-radius relations for rotating NSs. 
In both the cases, the maximum masses in presence of the magnetic field are 
less than that in absence of the magnetic field. 
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Figure 2: Mass-radius relations without rotation. The solid curve is for without magnetic 
field. The dashed, dotted and dot-dashed curves correspond to B c — 10 18 , 2 x 10 18 and 
3 x 10 18 G respectively. 



3. 3. Matter instability due to magnetic field 

It has been noted that there is an onset of matter instability for certain 
sets of and 7 (which correspond to different magnetic field profiles) for every 
B c . Now we explore EoSs for different magnetic field profiles to understand 
the instability. Figures [5^ andE^i show the variations of P at B c = 10 18 and 
3 x 10 18 G respectively, for different sets of j3 and 7. For each /3, we take 
minimum 7 to be 1 for which B varies very slowly with n& (see Figs. [5b and 
[6b)- The maximum value of 7 for every /3 has been chosen in such a way 
that beyond a certain value of n&, P ceases to increase (and decreases) with 
the further increase of rib- This implies that the matter becomes unstable 
after that value of n&. Keeping /3 constant if 7 decreases, then instability 
disappears at some value of 7 and the EoS tends to become that for without 
magnetic field. The point to note is that if (3 decreases, the instability occurs 
at larger values of 7 and rib, which is evident from Figs. Eh andEK- 

The instability occurs because the magnetic pressure contributes nega- 
tively to the matter pressure as reflected from the Eqs. (J21I) and (|25|) . Since 
the field strength increases with the increase of density, more negative contri- 
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Figure 3: Same as Fig. 2, except for rotating NSs with rotational frequency £1 = 2000 s 1 . 

bution is added to the pressure, and finally, at a certain density, the pressure 
ceases to increase with the increase of the density. This effect has not been 
discussed in any earlier literature to the best of our knowledge. This is con- 
sistent with the result of softening the EoS with the increase of magnetic 
field. Since the softening of the EoS implies that with the energy density the 
rate of increase of pressure decreases, above a certain value of magnetic field, 
the rate becomes zero, and subsequently becomes negative. That is why, 
above a certain value of magnetic field the matter becomes unstable. This 
constrains the value of central magnetic field of a NS having high surface 
magnetic field, e.g. magnetar. In this context, it is worth mentioning that 
the inclusion of AMM effect 3l| stiffens the EoS, which dominates over the 
effect due to Landau quantization on EoS when field strength is > 5 x 10 18 
G. This stiffening is only possible when the field pressure is not considered 
in obtaining the EoS. However, the inclusion of field pressure further softens 
the EoS, as discussed above. 

When the field strength is ~ 5 x 10 19 G, the deviation in EoSs between 
the cases with and without AMM effect is significant throughout the density 
However, this deviation is much small compared to that due to 
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range 

the extra effect arised from field pressure. Hence, at such a strong magnetic 



13 




Figure 4: Same as Fig. 2, except for rotating NSs with rotational frequency il = 4000 s 1 . 

field, AMM effect is unable to stiffen the EoS. It is, however, seen that when 
the field strength is near 5 x 10 18 G, AMM effect is insignificant at high 
density, but significant at low density. However, one should also remember 
that this result is obtained by considering constant field strength over the 
whole range of density. Thus, in the realistic case, when the field strength is 
several orders lower than 10 18 G at low density, i.e., near the surface, AMM 
effect does not play any role in modifying the EoS. Therefore, we do not 
consider the AMM effect in our calculation while B c < 10 19 G. 

In Fig. the mass-radius relations for different B c as well as that without 
magnetic field have been shown. We notice that for a fixed magnetic field 
profile, e.g. (3 = 0.1 and 7 = 1, the maximum mass decreases as B c increases 
from 10 18 G to 4 x 10 18 G. Now as B c increases further from 4 x 10 18 G, with 
the same magnetic field profile (same fl and 7), P starts decreasing with the 
increase of rib at a very low rib, leading to unstable matter. For B c = 5 x 10 18 
G, we obtain stable matter EoS for /3 = 0.01 and 7 = 2, which deviates 
maximum from that without magnetic field, as shown in Fig. |HJ It is found 
that the maximum deviation of mass-radius relation for B c = 5 x 10 18 G with 
respect to that without magnetic field is less than that deviation for B c = 
4x 10 18 G (see Fig. [7]). Now if B c increases further to 10 19 G, we note that (see 
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Figure 5: Variation of (a) P, and (b) B with normalized baryon number density for 
different magnetic field profiles with B c = 1 18 G. The solid curve is for without magnetic 
held. The sets of dashed, dotted, dot-dashed and short-dashed curves are for /3 — 0.1, 
0.01, 0.001 and 0.0001 respectively. In each set (for each (3), the upper curve in (a) is for 
7 = 1. The values of 7 for the lower curve in sets in (a) are given in Table [21 In (b), for 
each (j, the curves which reach B c at a comparatively lower n& correspond to maximum 7 
cases. 
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Table 2: 7-s for different lower curves shown in Figs. EK andEK for different /3-s. 



Fig. EJ) with the above said profile (ft = 0.01 and 7 = 2) the matter becomes 
unstable at a very low density. We have stable matter up to a reasonably 
high density with B c = 10 19 G for (3 = 0.001 and 7 = 3 with maximum 
deviation in EoS from that without magnetic field. However, Fig. [7] shows 
that in this case the deviation of EoS from that without magnetic field case 
is even lesser. Hence we can conclude that the effect of magnetic field is 
significant when B c is around 4 x 10 18 G depending on the magnetic field 
profile. Beyond that, the magnetic field profiles needed to create significant 
effect of magnetic field on the mass-radius relation make the matter unstable. 
As we have mentioned earlier, the effect of magnetic field on quark matter 
using the MIT bag model has been studied in many earlier work. Therefore, 
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Figure 6: Same as Fig. [5] except for B c = 3 x 10 18 G. 



further it will be interesting to see the effect of magnetic field on other models, 



e.g. of Dey et al. 29], which is in progress. 



4. Summary 

We have studied the effect of high magnetic field on the NS matter. The 
nuclear matter inside the NS is highly dense (2 — 15 times normal matter 
density). We have considered the possibility of appearance of hyperons at 
the core of the NS with the increase of density. Our model is based on the 
relativistic mean field theory for nuclear matter. With the increase of density, 
the heavier particles appear gradually. However, due to repulsive potential, 
S-hyperons do not appear. 

Observationally, the inferred surface magnetic field of a NS may be as 
high as 10 15 G. It is believed that at the center, the magnetic field is higher 
than its surface value. As the density of matter also decreases with the 
increase of radius, we have taken the parametrization of the magnetic field 



as a function of density J49J. We have studied the effect of magnetic field 
on the NS matter with the different sets of parameters j3 and 7. We have 
found that the effect of magnetic field is important when B c > 10 17 G. On 
the other hand, when B c ~ 10 19 G, the matter becomes unstable. This gives 
a maximum bound of B c for NSs. For most of the X-ray and radio pulsars, 
the surface magnetic field ranges from 10 8 to 10 12 G. It can be assumed that 
their central magnetic field is larger than their surface values by 3 — 4 orders 
of magnitude. Therefore, for this class of objects the effect of magnetic 
field is not significant. The effect of magnetic field can be observed only 
for magnetars whose surface magnetic field is 10 14 — 10 15 G and the central 
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Figure 7: Mass-radius relations for static star with different values of B c . The solid curve 
is for without magnetic field. The dashed curves are for B c = 10 18 , 3 x 10 18 and 4 x 10 18 
G respectively from top to bottom for j3 — 0.1 and 7=1. The dotted curve corresponds 
to B c = 5 x 10 18 G for (1 = 0.01 and 7 = 2, and the dot-dashed curve is for B c = 10 19 G 
when /3 = 0.001 and 7 = 3. 



magnetic field is constrained by the maximum value ~ 10 19 G. In reality 
also, we do not observe any compact object with ultra-strong magnetic field. 
This is consistent with the observed surface value of magnetars. 

It should be noted that in strong magnetic field, there may arise anisotropic 
effect due to the magnetization of matter. However, in NS matter the mag- 
netization is small [21]. In none of the previous work, this effect has been 
included. Therefore, it is expected that the inclusion of anisotropic effect 
might alter the present result quantitatively, but not qualitatively, leaving 
the conclusions unaltered. However, we have considered the field pressure 
contribution to the EoS of NS matter, which was not considered in any pre- 
vious work. 

In the presence of strong magnetic field, the nature of interaction between 
baryons may or may not be affected. However, the effect of magnetic field 
on nuclear interaction at very high density is poorly known till date, while 
nuclear interaction itself at high density is not well understood. Hence, we 
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Figure 8: P as a function of normalized baryon number density with B c — 5 x 10 18 G for 
different magnetic field profiles. The solid curve is for without magnetic field, the dashed 
curve for (3 = O.f and 7 = f, the dotted curve for j3 = O.Of and 7 = 3, the dot-dashed 
curve for f3 = 0.01 and 7 = 2 and the short-dashed curve for j3 — 0.001 and 7 = 4. 



have employed the same model for baryonic interaction, what is widely used 
for nuclear matter without magnetic field. As our main aim is to study the 
properties of NS matter in presence of magnetic field, we do not study the 
magnetic field effect on the NS structure. It appears that the NS/magnetar 
structure with realistic EoS in presence of the magnetic field depends on the 
radial distribution of the magnetic field. Hence, a more realistic treatment 
including the effect on NS structure due to anisotropy arised from general 
relativity is expected to strengthen the present conclusions. 



Acknowledgment s 

This work is partly supported by a project, Grant No. SR/S2HEP12/2007, 
funded by DST, India. The authors would like to thank Ritam Mallick of 
IISc for verifying one part of the computations. 



140 



> 




n b /n 



Figure 9: P as a function of normalized baryon number density with B c = 10 19 G for 
different magnetic field profiles. The solid curve is for without magnetic field. The dashed 
curve is for j3 = O.Of and 7 = 2, the dotted curve for f3 — O.Of and 7 = 1 and the 
dot-dashed curve for f3 — O.OOf and 7 = 3. 
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